A lattice model for the SU(iV) Neel-VBS quantum phase transition at large iV 
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We generalize the SU(TV = 2) S — 1/2 square-lattice quantum magnet with nearest-neighbor 
antiferromagnetic coupling (Ji) and next-nearest-neighbor ferromagnetic coupling (J2) to arbitrary 
TV. For all TV > 4, the ground state has valence-bond-solid (VBS) order for J2 = and Neel 
order for J2/J1 2> 1, allowing us access to the transition between these types of states for large TV. 
Using quantum Monte Carlo simulations, we show that both order parameters vanish at a single 
quantum-critical point, whose universal exponents for large enough TV (here up to TV = 12) approach 
the values obtained in a 1/TV expansion of the non-compact CP^ -1 field theory. These results lend 
strong support to the deconfined quantum-criticality theory of the Neel-VBS transition. 



Just as the destruction of magnetic order by thermal 
fluctuations is the paradigmatic example of a classical 
critical point pQ, the destruction of magnetic order at 
T = by quantum fluctuations is a prototypical example 
of a quantum-critical point [2 . Sometimes the quan- 
tum case can be entirely different due to novel quan- 
tum interference effects, which have no natural classical 
analogues. An important example is that of the square- 
lattice SU(TV) quantum antiferromagnet, where the de- 
struction of SU(TV) symmetry-breaking Neel order in the 
background of uncompensated Berry phases results in a 
valence-bond-solid (VBS) state with broken translational 
symmetry [3 - 5] . Recent work has presented speculative 
albeit compelling arguments that a direct generically con- 
tinuous Neel-VBS transition can exist. At such a decon- 
fined quantum critical (DQC) point both order param- 
eters are simultaneously critical [6 , a striking feature 
which is not contained in the conventional field-theory 
description of two independent order parameters (where 
one would instead expect a generic direct transition to 
be first-order). 

Given the major paradigm shift that could be spawned 
by the DQC idea, it has been of great interest to verify 
its validity by unbiased numerical studies of lattice spin 
models (Hamiltonians) that harbor Neel-VBS transitions 
[7KTT] . The weight of evidence from such work indicates 
that generically continuous Neel-VBS transitions indeed 
exist in quantum spin systems, with initial skepticism 
HI] appearing increasingly unfounded [13] [H] ■ 

This DQC scenario predicts that the SU(TV) Neel-VBS 
quantum-critical point falls into the universality class of 
the (2 + l)-dimensional non-compact CP^ 1 field theory 
[Pol [Ttlj . The connection between the phase transition in 
a microscopic Hamiltonian and the low-energy continuum 
theory description relies on speculative assumptions that 
are yet to be demonstrated convincingly. In order to 
provide support for this connection, therefore, one has 
to compare universal properties arising from these two 
starting points. Currently, the only technique with which 
the properties following from the CP^ -1 action can be 
studied analytically is the 1/TV expansion [17], but it is 
not clear whether the results of this approach are valid 




FIG. 1. (Color online) (a) Black (white) lattice sites indicate 
the A (B) sublattice on which spins of the J1-J2 model |T]) 
transform as the fundamental (conjugate) representation of 
SU(TV). Ji couples nearest neighbors with an SU(TV) singlet 
projection and Ji connects next nearest neighbors with an 
SU(TV) permutation, (b) Phase diagram of the model as a 
function of TV and g = J2/J1. 

down to the most interesting case of TV = 2. To test the 
DQC theory in an unbiased manner, it is thus of utmost 
importance to find quantum models in which the SU(TV) 
Neel-VBS transition can be studied for arbitrary large 
TV, to compare the critical exponents with those of the 
1 /TV expansion of the continuum field theory. 

Until now, the Neel-VBS transition could be accessed 
in an unbiased manner only for TV < 4, by quantum 
Monte Carlo (QMC) simulations of the so-called J-Q 
model [7HTTI El], in which the S = 1/2 Heisenberg 
(J) model is supplemented by certain multi-spin interac- 
tions (Q) favoring a VBS state. Is clear that this model 
alone cannot access the quantum transition for larger- 
TV, however, simply because for TV > 5 the J model 
is itself VBS ordered [TH [05] and the Q term only in- 
creases the strength of this order. To remedy this prob- 
lem, we here introduce an SU(TV) symmetric general- 
ization of the J\-Ji Heisenberg model model, with an- 
tiferromagnetic nearest-neighbor coupling J\ and ferro- 
magnetic next-nearest-neighbor coupling J2 [illustrated 
in Fig. Qa)]. For all TV > 5, this model harbors a VBS 
phase for J2/J1 = and a Neel phase for J2/J1 = 00. 
By detailed unbiased QMC studies for 5 < TV < 12, we 
find that that the two phases are separated by a single 
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phase transition, with no signs of discontinuities even on 
the largest systems sizes studied (L x L spins with L 
up to 128). Most remarkably, the anomalous dimensions 
of the Neel and VBS correlation functions of the model 
for large- -/V shows quantitative agreement with the ana- 
lytically known [2"D1 - |2"2"] scaling dimensions from the 1 /N 
expansions of the non-compact CP^ -1 model. 

The J\-Ji model. — Our SU(7V) symmetric model is de- 
fined with a local Hubert space of N states on each site of 
the square lattice illustrated in Fig. [l|a). We adopt the 
representation used previously in both analytic [1] and 
numerical [HI works on bipartite lattices, where the 
sublattice-A states transform under rotations with the 
fundamental representation of SU(iV), and the B sub- 
lattice states transform with the conjugate of this rep- 
resentation; \a) A -> U ap \P) A , \oc)b -> U* p \i3)b. The 
state \a) A \a)B is, thus, an SU(JV) singlet. is de- 
fined to be the projector onto this singlet between two 
sites i and j on different sublattices, i.e., H+j = —Pij/N 
is the SU(iV) generalization of the familiar Heisenberg 
antifcrromagnetic exchange (up to a constant). An- 
other simple SU(iV) invariant interaction is the permuta- 
tion operator between two sites on the same sublattice, 



Tlij\af3) = \/3a), so that Hi 



-Hij/N is the generaliza- 



tion of the the familiar ferromagnetic Heisenberg inter- 
action. The Hamiltonian we study here is given by 



H 



(1) 



where (ij) and ((ij)} denote first (A-B) and second (A-A 
and B-B) neighbor sites, respectively. 

With J 2 = it is now well known that the J\ model 
is Neel ordered for N — 2,3,4 and develops VBS order 
for N > 5 [151 AS ] - On the other hand, with Ji = each 
sublattice forms a trivial ferromagnet. A small J\ <C J% 
will clearly lock the individual sublattice magnetizations 
into a collective Neel ordered state. Thus, for each N > 5 
there must be an intermediate value of g = J2/J1 at 
which there is a quantum transition between these two 
phases (as we do not expect any other intervening phase) . 

QMC simulations. — All off-diagonal matrix elements 
in Eq. ([ij are explicitly negative and, hence, the model 
is free of QMC sign problems [and it also satisfies Mar- 
shall's sign criterion, ensuring an SU(iV) singlet ground 
state] . To obtain exact (within statistical errors) numer- 
ical results for its properties on large L x L lattices, we 
use the stochastic series expansion QMC method with 
global loop updates [2"3Tl25| . Throughout this work, we 
set J\ — 1 and the inverse temperature /3 — L/Ji (re- 
flecting the expected [5] dynamic exponent z = 1). 

We characterize the Neel phase as one with a finite spin 
stiffness p s (measured by the fluctuations of the winding 
number W of world lines; f3p s = (W 2 ) [351 023)- I n the 
magnetic phase, the static (oj = 0) Neel order-parameter 
susceptibility xn diverges with the "quantum volume" 
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FIG. 2. (Color online) Curve crossings used to locate the 
critical point for magnetic [VBS] order in the SU(5) J1-J2 
model are shown in the left [right] panel. The quantity f3p s 
[£vbs/L] diverges in the magnetic [VBS] phase and goes to 
zero in the non-magnetic [non-VBS] phase when (3 — L/J\. 
At a point where magnetic [VBS] fluctuations are critical, fip a 
[fvBs/£] becomes L-independent. These properties result in 
crossings of curves for different L at the critical point. The 
width of the vertical line shows the range of estimates of the 
common Neel- VBS critical point; g c = 1.615(10). Fig.[3]shows 
the analysis of the crossing points giving this result. 

of the system according to xn ~ jil? . We define the 
SU(iV) VBS correlation function using the operator P de- 
fined above in Eq. ([!]); Cv(r,r) = (-P o .o+x(0)P r: r+x(T))- 
(Po,o+x(0)) 2 . When Fourier transformed at uj — 0, q = 
(-7T, 0) it gives xv ■ This quantity can be used to test 
for VBS order since it diverges in the VBS phase as 
Xv ~ fiL 2 . We also use the standard definition of the 
correlation length of the VBS order £y as the square root 
of the second moment of the spatial correlation function 
Cy • Using these quantities we tested for long-range Neel 
and VBS order as the ratio g — J2/J1 is varied for each 
N and arrived at the phase diagram shown in Fig. |ljb) . 
We elaborate on the quantitative analysis below. 

Nature of the phase transition. — Fig. [2] shows QMC 
results for j3p s and ^vbs/-^ as functions of the cou- 
pling ratio g for the SU(5) model on lattices of size 
L = 8, 16, 32, 64, and 128. The quantum-critical point for 
the magnetic and VBS orders can be located by analysing 
crossing points versus g in j3p s and £vbs/£, respectively, 
computed on two different system sizes. As is clearly 
evidenced directly from this raw data, there are cross- 
ing points within a narrow window of g for both Neel 
and VBS orders, and these crossing points drift toward 
a common value g c with increasing L. 

In Fig. [3] we have plotted the crossing points between 
L and L/2 curves of j3p s and £vbs/£ for SU(7V) sys- 
tems with V = 5, 6, 10, and 12. Numerical extrapola- 
tions of the crossing data for both Neel and VBS orders 
in the SU(5) and SU(6) cases (top two panels) provide 
compelling evidence that in the thermodynamic limit 
the crossing points for both order-parameters approach 
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FIG. 3. (Color online) Convergence to a common critical 
point of the finite-size estimates for the Neel and VBS phase 
transitions for SU(5), SU(6), SU(10), and SU(12). The cross- 
ing points of L and L/2 curves for the same quantities as in 
Fig. [5] are shown for Neel (blue circles) and VBS (red squares) 
order for L < 128. For all four cases, the data are consistent, 
within error bars, with both order parameters becoming criti- 
cal at the same point when L — > 00. Numerical extrapolations 
(when reliable) are shown by dashed lines. For SU(10) and 
SU(12) the VBS correlation function has significant sublead- 
ing corrections that can be ignored only for the largest sizes 
(see [27] ). The grey lines show our estimates of the critical 
points (with their widths corresponding to the error bars). 

a common critical point. For SU(10) and SU(12) (bot- 
tom two panels) there is non-monotonicity in the VBS 
crossing points, making reliable numerical extrapolations 
difficult. The source of this behavior is subleading cor- 
rections in the VBS correlation function that increasingly 
dominate the rapidly decaying leading power-law behav- 
ior as N increases (see note However, even for the 
SU(10) and SU(12) cases, for the largest system sizes 
the leading behavior still dominates the subleading cor- 
rections, and the VBS crossings are consistent with the 
same critical point as obtained from the better-behaved 
f3p s crossing in the L — > 00 limit. There is, thus, good 
reason to believe that both orders go critical at the same 
point for all N. We have found no evidence for hysteresis 
or double peaked histograms that would be expected for 
a first-order transition. 

Since a reliable numerical extrapolation of the /3p s 
crossing is possible for each N, we use this quantity to 
extract the quantum-critical points. Note that possible 
weak corrections to standard scaling behavior in quan- 
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FIG. 4. (Color online) Examples, for the SU(5) and SU(8) 
models, of power-law scaling of the Neel and VBS order pa- 
rameters close to the quantum-critical point. According to 
Eq. (|2| and standard finite-size scaling arguments, \n / ' f)L 2 ~ 
and Xv /f3L 2 ~ 1/L 1+Vv . The linear regression fits 
(dashed lines) allow us to extract nv and vn, which are shown 
in Fig. [5] as a function of N. 

tities whose scaling from depends only on the dynamic 
exponent, such as p s , which have been discussed for the 
J-Q model [TTJ [T3l [H] , would not affect the above anal- 
ysis of crossing points versus g (see Appendix D in |11|). 
Corrections to scaling of the quantity j3p s itself will be 
studied in detail elsewhere. 

At a common quantum-critical point, it is expected 
that the otherwise independent Neel and VBS order pa- 
rameters will both have correlation functions that decay 
as power laws in space (r) and imaginary time (r); 

C Ny (r, r) ~ (r 2 + c 2 r 2 )-( 1+ ^«)/ 2 . (2) 

The exponents 77 n and r)v, the so-called "anomalous di- 
mension" of the Neel and VBS fields, are universal num- 
bers according to standard theory of critical phenomena. 
Universality implies that they are independent of details 
of the microscopic interactions of the model from which 
the correlation functions are extracted, but they do de- 
pend on the symmetry of the model, i.e., in our case they 
should only depend on of the SU(A r ) symmetry. To es- 
timate these exponents with QMC simulations, we have 
studied the A^ = 5, 6, 8, 10, and 12 models at values of the 
coupling ratio g within the estimated critical points from 
the analysis shown in Fig. [3] We extract the exponents 
from the size dependence of the correlation functions, as 
illustrated for two cases in Fig. [4] 

Here we note again that the scaling behavior in the J-Q 
model of quantities depending only on the dynamic expo- 
nent show weak corrections to scaling [TTJ [T^l [H] . If such 
corrections are present also in the correlation functions 
analyzed here, they would somewhat affect the values of 
the exponents, e.g., a multiplicative log is hard to distin- 
guish from a small change in an exponent (see |25| for an 
example in one-dimension). 
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FIG. 5. (Color online) Anomalous dimensions of the Neel 
(left) and VBS (right) fields extracted from the critical scaling 
analysis. The main panels show vn and nv versus 1/JV. For 
N = 2, 3 and 4, the data are for the J-Q model [TO], and the 
results for JV > 4 are for our J1-J2 model. The analytic results 
from the 1/N expansion of the CP^ -1 field theory are shown 
as thick red lines. The left and right insets show JV(1 — n jv) 
and (l + 7)v)/N, respectively. These quantities must be finite 
in the JV — > 00 limit according to the DQC theory and should 
be given by Eq. |3| (solid straight lines in the insets). The 
next corrections to the exponents have not been computed 
analytically yet, but we can estimate them approximately as 
1 + Vv = 0.2492JV + 0.68(4), -q N = 1 + 32/(tt 2 JV) - 3.6(5)/JV 2 
(shown as dashed lines). 



Comparison with large-N results. — The extracted ex- 
ponents rjM and r\y are shown versus 1/N in the main 
panels of Fig.JH} For SU(10) and SU(12), r\ v becomes too 
large to extract reliably (see note [57] ) . Our main objec- 
tive is to compare the results with those of the CP^ -1 
universality predicted by the DQC scenario (§]. Analytic 
large- JV results currently available for these indices are 



r) N = 1 - 32/(tt 2 JV), l + nv = 2<5iJV, 



(3) 



where Si ~ 0.1246. It is interesting to note that the 
calculations underlying these two results are entirely dif- 
ferent. While rjN was obtained from an 1/JV expansion 
of the Neel order parameter expressed in terms of the 
CP^ 1 fields [20], rjv was computed [HJ[22] by exploit- 
ing a non-trivial relation, predating the DQC theory, be- 
tween monopoles in the field theory and the VBS order 
on the lattice @]. 

In Fig. [5] we show the results of the 1 /JV expansion, 
Eq. ([3]), as continuous curves. In the insets we plot the 
same data in such a way that in the 1/JV — > limit we 
can do a more direct comparison with the two irrational 
numbers 32/-7T 2 and 2<Ji that are predicted based on the 
Q-pN-i theory We have also fitted these data to straight 
lines, which give numerical predictions of the next-higher 
corrections to the large- JV forms in Eq. pi) . 

The most important features in Fig. [STthat lend sup- 
port to the DQC scenario are: (1) The exponents r]N 
and r\v for the JV > 4 J\-J% models are consistent in 



trend with previous estimates for JV = 2, 3, 4 based on 
the entirely different J-Q model [10 . This is line with 
the concept of universality between different microscopic 
models, characteristic of a continuous transition. (2) The 
Neel-order exponent ijn connects to the leading 1/JV be- 
havior, approaching 1 as JV — > 00. This large value (in 
contrast to the normal mean-field value 77 = and typi- 
cally very small values at conventional critical points) was 
one of the important early "smoking gun" predictions of 
the DQC scenario [B]. (3) The leading 1/JV correction to 
t)n is of the correct sign and within a few percent in mag- 
nitude of the analytic result 32/n 2 for the CP^ -1 theory 
[20] . (4) Vv increases rapidly with JV and has a trend 
that is fully consistent with the very non-trivial large- JV 
result [2UH22]- 

Conclusions. — We have provided an SU(JV) symmetric 
sign problem free model that allows, for the first time, 
unbiased studies of the Neel- VBS transition for arbitrary 
JV > 5 on large lattices. Our model generalizes to SU(JV) 
the important J1-J2 model for SU(2) spins, which has a 
long history of studies by exact diagonalization and var- 
ious approximate methods (normally for antiferromag- 
netic J2, in which case it has a Neel- VBS transition also 
for JV = 2, which, however, is difficult to study reliably 
due to QMC sign problems) [25] , Our model opens new 
avenues for detailed connections between large- JV theory 
and unbiased numerical simulations. 

In this initial study, we found direct Neel- VBS tran- 
sitions for 5 < JV < 12, with no signs of discontinuities 
up to large system sizes (L < 128). When analyzed as 
continuous quantum critical points, we fond remarkable 
agreement for universal exponents with the 1/JV expan- 
sion of the CP^ -1 field theory. This quantitative com- 
parison lends very strong support to the DQC theory [5J. 
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